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Abstract 


A  distributed  parameter  model  for  a  class  of 
search  and  capture  problems  In  the  plane  Is  derived. 
The  output  Is  a  time-varying  simple  loop  In  R2  and  the 
Input  a  moving  line  segment.  The  applications  dis¬ 
cussed  Include  destroyer  search,  oil  spill  clearance, 
forest  fire  control  and  fishing.  Circular  motion  Is 
shown  to  produce  a  non-trlvlal  steady-state  region 
whenever  the  velocity  ratio  between  searcher  and 
evader  attains  a  computed  critical  value. 

Introduction 

Suppose  that  we  wish  to  carry  out  a  search  In  the 
plane,  R2,  with  the  visible  points  at  time  t  being  a 
line  segment  of  fixed  length  d  attached  to  a  moveable 
origin  0.  (See  Figure  l-a.)  The  search  might  begin 
with  a  rotation  of  "the  line  of  sight"  as  Indicated  In 
Figure  l-b. 


Fig.  I  Starting  the  search 

The  objective  is  to  search  for  evasive  objects  In  the 
vicinity  of  0  which  have  speed  known  to  be  bounded  by 
a  given  number  B  >  0.  If  after  a  steady  rotation  of 
the  line  of  sight  by  360°  no  evaders  have  been 
detected  then  we  can  be  assured  there  are  no  evaders 
present  In  the  clear  region  Indicated  In  Figure  2. 
Notice  that  the  clear  region  Is  a  proper  subset  of  the 
searched  region  since  the  objects  sought  can  move 
about. 

If  the  search  were  now  terminated  the  clear  region 
would  tend  to  contract.  To  counter  this  tendency  with 
the  hope  of  eventually  expanding  the  clear  region  the 
searcher  might  proceed  with  some  combined  motion  of  0 
and  the  line  of  sight.  The  partial  differential  equa¬ 
tion  derived  In  this  paper  provides  a  mathematical 
foundation  for  studying  the  relative  merits  of  such 
search  procedures. 

Other  game-theoretic  problems  to  which  the  model 
also  applies  are  discussed  In  the  later  sections.  A 
problem  In  which  the  searcher  moves  along  a  circular 
path  Is  studied  In  some  detail. 

Some  reflection  upon  the  nature  of  the  search 
situation  presenting  Itself  In  Figure  2  suggests  that 
the  appropriate  choice  for  the  state  of  the  system 
might  be  a  time-varying  loop  bounding  the  searched-and- 
clear  region.  The  control  Input  will  be  the  moving 


Fig.  2  The  searched-and-clear  region 

assuming  no  detection  has  occured 

straight  line  segment  extracted  from  the  line  of  sight 
by  Intersection  with  the  loop.  This  segment  will 
hereafter  be  called  the  control  or  Input  segment. 

Derivation  of  the  model 

To  each  simple  loop  of  the  type  eluded  to  we 
attach  coordinates  In  the  manner  Indicated  in  Figure  3. 


Fig.  3  Assignment  of  system  coordinates 

Looking  at  the  loop  at  time  t  first  locate  the  control 
segment  which,  when  extended,  penetrates  the  Interior. 
Along  this  extended  line  select  an  origin  0.  The 


3 


state  of  the  system,  l.e.  the  loop,  can  then  be 
described  by  a  function  p(t,$)  where  $  Is  the  Indi¬ 
cated  angular  coordinate.  To  describe  the  motion  of 
the  control  segment  a  control  point  such  as  the  center 
Is  selected  along  It.  The  radial  distance  to  the  end 
of  the  loop  nearest  0  Is  denoted  by  r(t)  and  fi(t)  Is 
the  segment's  angular  velocity.  The  respective  speeds 
of  the  control  point  along  the  segment  and  orthogonal 
to  It  are  Indicated  as  x  and  y.  At  any  time  that  a 
point  on  the  loop  Is  not  being  driven  by  the  control 
segment  It  will  be  assumed  to  be  moving  Inward  along  a 
normal  to  the  loop  at  speed  8. 


Looking  ahead.  It  turns  out  to  be  convenient  to 
describe  the  system  In  terms  of  the  following  canon- 


u4  =  x 


variables. 

(1) 

p(t,O)0  -  r8 

(2) 

r  -  p+(t,0) 

(3) 

(4) 

r 

(5) 

and  to  position  the  control  point  at  the  outer  end  as 
Indicated  In  Figure  4.  There  u_  and  u,  are  identified 
as  the  respective  speeds  of  0  orthogonal  to  and 
parallel  to  the  moving  segment.  The  variables  u.  and 
u_  give  the  speeds  of  the  ends  of  the  control  segment 
along  the  line  of  sight. 


Observe  that  the  change  of  control  variables  (l)-(5) 
can  be  Inverted  by  Integrating  that  system  of  linear 
ordinary  differential  equations.  Upon  assigning 
rectangular  coordinates  (z.,z,)  to  the  origin  of  the 
loop  It  Is  clear  that 

Zj  ■  UjCose  -  u2slne  (6) 

z2  *  ursine  +  UjCOsB  <7) 

which  points  out  that  the  motion  of  0  Is  strictly 
governed  by  U|,  u2  and  Uj. 

In  Figure  5  we  study  the  change  In  the  geometry 
or  the  system  during  a  small  time  displacement  At.  To 
each  point  on  the  loop  at  time  t  with  coordinates 
♦  ,  o(t,4)  there  is  a  corresponding  point  with  coordi¬ 
nates  ♦,  p(t+At,f)  on  the  displaced  loop.  This 
correspondence  produces  the  Induced  angles  At,  y  and 
n  Indicated  In  the  diagram  along  with  the  magnitude 
of  the  loop  displacement  para  I  led  to  Itself,  0At. 


Fig.  5  The  change  In  geometry  during  At 

First  we  compute  4,  the  Induced  change  In  p 
Indicated  in  Figure  5.  Consider  the  Infinitesimal 
right  trlange  generated  by  the  displacements  4,  BAt 
and  replacement  of  the  loop  at  the  point  A,  p(t+At,A) 
by  its  tangent  line.  (See  Figure  6.) 


Fig.  6  The  Infinitesimal  triangle 
for  computing  4 

Recall  that  In  the  elementary  differential  geometry  of 
a  curve  with  polar  coordinates  0,  r (6 )  (not  to  be 
confused  with  tho  control  variables)  the  standard  angle 
A  between  the  radial  and  tangent  line  Is  given  by 

tent  ■  .  Thus  we  compute 

r6 

4  ■  Jim)2 1  i2(M)  ♦  p|(t,»).  (co 

v  tan>  p(t,A>  *  ♦ 

Next  we  apply  the  law  of  sines  In  Figure  5  to  got 


p(ttAt.A)  +  4  _  sln(n-A-AA) 
pTTTi+M)"'  slniyta) — 


(9) 


By  adding  and  subtracting  the  appropriate  terms  this 
last  equation  can  be  written  as 


6  sin(n~$~A$ 
XT  sinfy+t) 


>  A£  f  1‘JJjlI' h}2  ~ 

-  At  L  "  *  H  J 


I  [s  In  (n-i-A*)  .1  ..  .  rp(tfAt.*)-|.(f.»)1 

'  XT  [  jln(y+J )  'J  ^ u-j 

from  figure  5  It  is  clour  tti.it 


sin(n-t-A$) 

sln(y+i) 


+  I  as  At  +  0. 


Notice  that  y  -  A0  +  n  =  it.  Hence  y  *  it  ■>  n  as  it  ♦  0 
and  by  the  employment  of  standard  trigonometric 
identities 


I  T sin(n-$-A$)  _  _  J _ sln(o-t-At)  _ 

At  [  sln(y+i )  j  ’  At  [sin(n-^-A61  J 

2  /A9-Aii\  [C0S("'<I’~  A^A6)] 

=  At  sin( — i-^)  [sTn7rH>=^8T - J 

I  Im  (^2.  -  M)  r  |+tan(  I  im  n)tan  *1  ..  . 

*  Tan~TTTrrTr^»T  =  [ fan’{’ri¥" nl -"f an"‘rj 1  'mlI¥  ~  a£K 


But  tan ( I i m  n )  =  —  and 
u3 


im  |i=  a,  + 


u_cos$  *  u  *i  I  n<> 

/  J 


into  account.  For  zoro  ovader  speed  (B  =  0)  that  term 
drops  out  loavinq  a  guasi-l inoar  equation.  Its  role 
In  the  damping  of  tlio  solution  will  become  apparent. 


In  the  model  u. 


u^  are  distributed  controls  with 


u^,  Uj  strictly  boundary  controls  and  Uj  mixed.  The  u( 

variables  might  well  be  called  Input-output  variables 
(distinguished  from  the  state-the  loop)  since  from  the 
geometric  nature  of  the  model  It  is  clear  they  can  not 
be  chosen  arbitrarily.  The  Input-output  nature  of 
u^,  u5  depends  upon  whether  Uj,  u^,  Uj  are  In  the  mode 

tending  to  stretch  or  compress  the  ends  of  the  loop. 

The  model  will  exhibit  memory-loss  If  the  control 
segment  Is  "backed-up"  over  the  Interior  of  the  loop. 
This  Is  a  sacrifice  made  to  maintain  the  relative 
simplicity  of  the  model.  It  is  also  capable  of 
exhibiting  shock  and  bifurcation  phenomena. 

Game-theoretic  applications 

The  originally  stated  goal  underlying  the  deri¬ 
vation  of  ( 15)— < 18)  was  a  model  for  generating  the 
boundary  loop  of  the  searched-and-clear  set  for 
search  in  which  the  field  of  visibility  at  any  given 
time  Is  a  line  segment  and  the  objects  being  sought 
have  evasive  speed  bounded  by  a  parameter  B.  We  shall 
point  out  some  other  conflict  situations  to  which  the 
model  applies  as  well.  The  control  theory  of  this 
distributed  system  will  be  presented  In  another 
article.1 


both  direct  conclusions  drawn  from  Figure  5.  Substi¬ 
tution  Into  (12)  then  shows 


I  |  sin(n-$-A<f ) 
At  LsTnTy+jT 

as  At  -»  0. 


CujSlnif  t  UjCOS*] 


By  substituting  (8)  Into  (10)  and  then  using  (II),  (13) 
and  (14)  to  take  limits  as  At  -*  0  the  partial 
differential  equation  In  the  following  model  Is 
deri  ve.J. 

The  distributed  model 


t  =  [ul  +  ~ 


jcosfi  -  u.sinf 


-  u_sin$  -  u,cos$  -  -  Vpz  t 
z  j  p  ♦ 

p(0,t)  =  p  (♦) 

0 

t 

p(t,0)  =  Po(0)  +  /Cuj(o)  -  Uj(o)]do 


Destroyer  search 


Suppose  that  an  airplane  flying  about  an  aircraft 
carrier  has  two  observers  on  board  each  watching  out 
one  of  the  side  windows  for  the  purpose  of  detecting 
any  enemy  destroyers  before  the  destroyers  can  bring 
the  carrier  within  range  of  their  missiles.  If  we 
wish  we  can  assume  that  the  dostroyers  have  superior 
radar  or  perhaps  satellite  Information  so  they  would 
maneuver  In  an  attempt  to  approach  within  range  of 
the  carrier  while  remaining  outside  the  detection 
range  of  the  radar  on  the  carrier  and  sight  of  the 
airplane.  Using  the  search  path  selected  by  the  air¬ 
plane  as  Input  ( 1 5 >— ( 1 8 >  would  generate  the  time- 
varying  region  which  the  carrier  could  be  assured 
contains  no  destroyers  as  long  as  the  airplane  reports 
no  sightings.  Clearly  the  airplane  would  want  to  choose 
a  path  which  would  expand  this  region  Into  a  comfortable 
area  about  the  carrier  (say  larger  than  the  destroyer 
missile  range  if  possible). 

Here  we  will  consider  only  circular  search  and 
look  for  the  path  with  the  largest  steady  state  disk. 
That  Is  we  select  U|(t)  =  C|,  constant,  and  u((t) 

(i  *  2,3)  Identically  zero  In  (15).  This  gives  rise  to 
the  ordinary  differential  equation  for  the  steady  state 


p(t,2n )  =  oq (2n )  t  /Cu^to)  -  Uj(o)]do 


The  initial  function  p  (•)  In  (16)  Is  prescribed. 

The  last  two  boundary  conditions  arise  from  the  choice 
of  positioning  the  control  point  (x,y)  at  the  end  of 
the  control  segment;  l.e.  requiring  that 
r(t)  »  p (t ,2i)  -  p(t,0)  In  the  Integrated  form  of 
r  =  u4  -  u5  and  p^ft.O)  *  Uj  -  Uj  which  follow  from 

(3),  (4)  and  (5). 

Remarks  about  the  model 

The  partial  differential  equation  (15)  contains  a 
nonlinear  term  which  reflects  the  evasion  being  taken 


c|%  •  p  VP2  +  f>\  -  o 

which  can  be  Integrated  to  give 

..fJGW- 


Equation  (20)  determines  a  one-parameter  family  of  . 
monotone  loops  p  «  p(t)  with  paramter  p(0)  »  p„  >  *•  . 

(See  Figure  7.)  Define  A  *  p(2w)  -  pQ,  the  length 
of  the  control  segment,  and  let  the  airplane  move 


along  at  the  center  of  the  control  segment  with  linear 
speed  v. 


c,  = 


°o  +  7 


(21) 


which  can  be  substituted  Into  (20)  along  with  $  *  2n 
■nd  p  *  <>o  ♦  A  to  get 


Tha  following  theorem  states  that  (22)  Implicitly 

y 

defines  s  p  as  a  function  of  the  parameters 

o  o  p 

1  and  ^  . 

Theorem  I.  If  j  Is  sufficiently  large  so  that 

V(|^>2  -  I  t  sin”1  <|^)  >  (23) 


I  ).9n : 

(a)  for  each  A  >  0  (22)  has  a  unique  solution 
>  >  0.  Moreover  (23)  Is  necessary  for  such  a 

.olutlon  to  exist. 


(b)  p  <A,-y)  ♦  •  as  4  +  •. 

(c)  po(a,j-)  ♦  0  as  4  ■*  0. 

Proof.  One  can  easily  verify  that  (23)  Implies 

that  ••  Hence  the  Integral  In  (22)  converges  to 

zero  as  p  ♦  *•.  On  the  other  hand  If  we  set 
0 

p  =  — -- —  In  the  Integral  and  apply  (23)  then 

we  see  the  Integral  takes  on  a  value  at  least  as  great 
as  2w.  The  existence  of  the  required  solution  p 

to  (22)  then  follows  from  the  Intermediate  value 
theorem.2  The  uniqueness  follows  from  the  fact  that 
the  Integral  Is  monotone  decreasing  In  po<  The 

necessity  of  (23)  Is  clear  from  the  monotonlclty  and 
the  observation  that  If  p  were  chosen  any  smaller 
than  the  above  choice  theR  the  Integral  would  not  be 
defined  since  the  lower  limit  of  Integration  would  be 
less  than  one. 

To  prove  (b)  note  that  If  for  some  sequence 
A.  •*  *  the  corresponding  values  of  p  remained  boundiM 
tnen  (22)  would  be  violated  since  thR  lower  limit  of 
Integration  would  again  become  less  than  one.  Part 
(c)  follows  from  the  observation  that  If  for  some 
sequence  A.  converging  to  zero  the  corresponding 
values  of  p  remained  bounded  away  from  zero  then  the 
limits  of  iRtegratlon  In  (22)  would  both  converge  to 

and  that  equation  would  be  violated. 

Remark.  In  the  destroyer  search  problem  j 

Is  the  observers'  range  of  visibility.  Theorem  I 
Implies  that  If  the  search  plane  can  attain  or  exceed 
the  critical  speed  ratio  (abocr  3.9)  then  once  the 
disk  of  radius  p  computed  from  (22)  has  been  cleared 
about  the  carrleP  the  carrier  enn  rest  assured  that 
there  are  no  destroyers  within  distance  PQ  from  It  as 

long  as  the  search  plane  reports  no  sightings.  More¬ 
over  this  Is  the  largest  such  disk  for  the  given 
search  path. 

Oil  spill  clearance 

Many  seemingly  different  conflict  situations  have 
essentially  the  same  underlying  mathematics.  For 
example  the  model  ( 1 5 >— ( 1 8 )  would  apply  as  well  to  the 
problem  of  cleaning  up  an  oil  spill  on  the  surface  of 
the  ocean  as  It  tends  to  dissipate  In  all  directions 
at  a  speed  6  or  less  due  to  wave  action,  wind,  etc. 
Here  we  are  thinking  of  a  clean-up  device  which  absorbs 
a  swath  of  the  oil  as  It  Is  steered  around  the  peri¬ 
meter  of  the  spill.  In  this  problem  6  would  be  a 
negative  parameter.  The  results  given  for  the  circular 
search  problem  would  apply  directly  to  give  the  radlu 
of  the  largest  circular  spill  that  could  be  contained 
and  eventually  cleaned  up  by  a  piece  of  equipment  of 
speed  v  and  swath  width  A. 

Forest  fire  control 


The  distributed  model  Is  well  suited  for  analysing 
a  variety  of  forest  fire  problems.  For  example  con¬ 
sider  the  situation  wherein  en  airplane  Is  dropping 
fire  fighting  chemicals  about  a  town  to  protect  It 
from  a  surrounding  fire.  If  the  fire  had  not  yet 
contacted  the  treated  region  we  would  set  0  *  0  In 
(1 5)- (18)  to  predict  the  extent  of  the  protection 
which  could  be  developed  by  any  available  equipment. 

The  control  segment  In  this  problem  would  be  the  cross 
sectional  line  of  the  chemical  trail  as  It  struck  the 
surface.  Once  the  fire  reached  the  treated  region  a 


positive  value  of  8  would  be  used  to  activate  the  non¬ 
linear  term  accounting  for  the  spread  of  fire  across 
the  treated  area.  The  spread  of  an  Initially  localized 
fire  would  be  modeled  using  8  negative. 

Netting  fish 

Suppose  that  a  fishing  vessel  towing  a  net  sights 
a  large  school  of  fish  swimming  near  the  surface  but 
the  school  quickly  submerges  from  sight  upon  being 
pursued.  If  a  bound  upon  the  swimming  speed  of  a  fish 
Is  known  can  the  fishing  boat  catch  all  the  fist  In 
Its  net?  This  Is  another  Interesting  question  which 
can  be  asked  of  our  model.  Of  course  we  are  making 
the  assumption  that  the  fish  remain  at  a  sufficiently 
shallow  depth  so  that  we  can  keep  the  problem  two 
dimensional.  The  answer  will  depend  upon  the  speed  of 
the  vessel,  the  size  of  Its  net  and  the  Initial 
distribution  of  fish. 


time  t  -  I  the  Initial  loop  has  been  steered  down  to 
the  semi-circular  loop  indicated  In  Figure  8.  That 
loop  can  be  further  shrunk  down  to  a  point  according 
to  the  solution 


p(t,t) 


2(2  -  t)  cost 


0  it  < 


.  <_  t  <  2" 


during  I  i  t  <_  2  by  switching  the  controls  to  Uj(t)  = 

Ug(t)  =  I,  u.(t)  *  -I  on  I  i  t  i  2.  Hence  the  model 
continues  to  hold  In  spite  of  the  fact  that  the  origin 
of  the  loop  strikes  the  loop  Itself  and  then  sticks 
there. 


Ac  know  I edqement 


Consider  the  model  with  8  =  I,  u,(t)  =  = 

Uj(t)  a  0,  u4(t)  *  u5(t)  =  -I  and  Initial  function 
f  I  ,  0  i  t  i  " 


I.  COS  t  +  V4 


sln2t  ,  "  1  *  I  2" 
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Fig.  8  The  Initial  loop  and  the 
state  at  switching 

The  partial  differential  equation  then  Is 


v 


Its  solution  could  be  obtained  using  the  method  of 
characteristics  although  It  is  much  easier  to  obtain 
the  following  solution  directly  from  the  Initial  loop 
and  the  geometric  nature  of  the  problem. 


